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Abstract 

Calculations of various corrections to the g factor of Li-like ions are presented, which result in a sig- 
nificant improvement of the theoretical accuracy in the region Z = 6-92. The configuration-interaction 
Dirac-Fock method is employed for the evaluation of the interelectronic-interaction correction of order 
and higher. This correction is combined with the 1/Z interelectronic-interaction term derived within 
a rigorous QED approach. The one -electron QED corrections of first in a are calculated to all orders in 
the parameter aZ. The screening of QED corrections is taken into account to the leading orders in aZ and 
1/Z. 

PACS numbers: 12.20.Ds, 31.30.Jv, 31.30.Gs 
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I. INTRODUCTION 

Recent high-precision measurements of g factors of H-like carbon 111 and oxygen [Q] have 
provided a possibility for studying binding QED effects in an external magnetic field in these 
systems. The experiments were performed on single H-like ions confined in a Penning ion trap at 
low temperatures. A relative accuracy of 5 x 10^^° was obtained in Ref. I2I for the ratio of the 
electronic Larmor precession frequency ujl and the ion cyclotron frequency uJc, which is connected 
with the (yf -factor value by 

_ 5f |e| miop 
uOc 2 q m 

where e is the elementary charge unit, q is the charge of the ion, mion is the ion mass, and m 
is the electron mass. The experimental results are shown to be sensitive to one- and two-loop 
binding QED effects and to the nuclear-recoil corrections. Further progress is anticipated from 
the experimental side, as well as an extension of measurements to the ions with more than one 
electron. 

New perspectives for testing QED effects in g factors of highly charged ions motivated numer- 
ous theoretical investigations on this subject during the last years. We mention here numerical 

nnnnn 

calculations of the one-loop self-energy correction to all orders in aZ ||4l la 13. ISD, numerical 
||5|,|6|] and analytical lull evaluations of the one-loop vacuum-polarization contribution, an- 

alytical results for the nuclear-size correction flll^, and calculations of the nuclear-recoil effects 
OQQIl^ and two-loop binding QED corrections OH [3]. As a resuk of these studies, 
the present theoretical accuracy of the (7-factor values is several times better than that of the exper- 
imental results. An important consequence of the investigations of the g factor is the possibility 
to determine the electron mass from Eq. ([T]) by combining the theoretical (7-factor value and the 
experimental result for ujl/uJc. We note that the 1998 CODATA value for the electron mass OOll 
has an error of 2 x 10~^, which is 4 times larger than the experimental uncertainty of the result 
for carbon ^ and 3 times l^er than the one for oxygen A new determination of the electron 
mass presented in Refs. [Q, I3 Izit provided an improvement of the accuracy of the electron 
mass by a factor of 4. As a result, the 2002 CODATA value for the electron mass |Q] is derived 
mainly from the g factor of H-like ions. It is expected that in the future an extension of experi- 
mental investigations towards higher- Z ions could lead also to an independent determination of 

nfi 

the fine structure constant a HI 8112311 . 

The accuracy of theoretical values for the g factor of high- Z H-like ions is presently limited by 
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nuclear effects lisll^l- The uncertainty introduced by them in the high-Z region is comparable 
with the binding QED correction of second order in a. Since the nuclear effects do not allow an 
accurate theoretical description at present, this puts a serious obstacle on the way to improvement 
of theoretical predictions and to an identification of two-loop QED corrections in future experi- 



ments. However, it was recently shown 1125 . 



2611 that the uncertainty due to the nuclear effects can 
be significantly reduced in a specific combination of the g factors of H- and Li-like ions with the 
same nucleus, 

g' = 9{is)^2s~^gis, (2) 

where gis and g^is)^ 2s are the g factors of H- and Li-like ions, respectively, and the parameter ^ is 



calculated numerically, as explained in Ref. ^25"]. Numerical calculations show that the influence 
of the nuclear effects on the difference g' is by two orders of magnitude smaller than that on the 
5f-factor values gis and g(isY 2s separately. Therefore, the specific difference g' can be in principle 
studied up to much higher accuracy than the g factor of an H-like or Li-like ion. In order to realize 
this project, one need to perform theoretical and experimental investigations of the g factor of 
Li-like ions with an accuracy comparable to that for H-like ions. 

Extending theoretical description from an H-like to a Li-like ion, one encounters a serious 
complication due to the presence of additional electrons. Moreover, all contributions to the g — 2 
value for an ns state are of pure relativistic and QED origin, which makes the treatment of the 
electron correlation much more intricate than, e.g., in calculating the binding energies. A number 
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of relativistic calculations of the g factor of Li-like ions were carried out previously 1127 

However, in order to reach the accuracy comparable to the one for H-like ions, a 



M, we 



systematic treatment of the QED corrections is required. In our previous publications 1125 . 
presented theoretical values for the g factor of Li-like ions including all presently known QED, 
nuclear, and interelectronic-interaction corrections for a wide range of the nuclear-charge number 
Z. The goal of the present investigation is to improve these theoretical predictions by calculating 
several corrections that provide the leading uncertainty to the theoretical values. 



For low-Z ions, the uncertainty of the (^-factor values II25L 12611 was mainly determined by the 
error due to the interelectronic-interaction correction. A part of this correction that is of first order 
in 1/Z was evaluated in Ref. l25i within a rigorous QED approach, whereas the remainder (of 
order l/Z"^ and higher) was extracted from the variational calculations by Yan in case of 

low-Z ions and from the many -body perturbation theory calculations by Lindroth and Ynnerman 



1 29|] for high-Z ions. In the present work, we evaluate the interelectronic-interaction correction of 



order and higher by a large-scale configuration-interaction Dirac-Fock method. 



For middle- and high-Z ions, the main error of the theoretical values of Ref. 112511 originated 
from the binding QED correction. In that work, it was evaluated to the leading order in oiZ 
[namely, {aZy] and to the zeroth order in 1 /Z. In the present investigation, we employ the recent 
calculation of the self-energy term of the one-electron QED correction to all orders in jsl 
and calculate the vacuum-polarization term. We also evaluate the screened QED correction to its 
leading orders maZ and 1/ Z. 

The relativistic units {h = c= V) and the Heaviside charge unit {a = e^/ (47r), e < 0) are used 
in the paper. 

II. FORMULATION OF THE PROBLEM 



The g factor of an atom with a spinless nucleus is defined as 

1 dEiU) 



d 



(3) 

H=0 



n=o f^oJ 9T-C 

where H is the Hamiltonian of the atom in magnetic field J-CH = |?i | and the z axis is directed 
along Ti, E{H) is the energy of the state with the maximal (Mj = J) projection of the angular 
momentum on the z axis, and /io = |e|/(2m) is the Bohr magneton. If the perturbation theory of 
first order in H is used to obtain the energy shift AE{T-C), Eq. Q can be written as 

^ ^ 1 AEjH) 

In what follows, we will use both definitions. 

In the case of a high-Z Li-like ion, in the zeroth-order approximation one can neglect the 
interaction of the 2s electron with the closed (Is)^ shell. We thus use the one-electron Dirac 
equation as the starting point in our evaluation of the g factor. The corresponding Hamiltonian is 

h = cx-p + l3m + V + 1/"^^^^ , (5) 

where V = —aZ/ris the Coulomb potential induced by the nucleus and V""^^^ = —ect-A 
represents the interaction with the magnetic field. For the homogeneous magnetic field the vector 
potential is A = l/2[n x r] and thus V'^^sn = _(e/2)?i ■ [r x a]. 

Evaluating the energy shift as the expectation value of V""^^^ with the Dirac wave functions, 
we obtain the lowest-order (^-factor value for the 2 s state, 

g^ = -(l + ^2 + 27] = 2 - + . . . , (6) 
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where 7 = a/1 — {aZ^. Various corrections to gB arise due to the interelectronic interaction 
(A^fint), one-electron QED effects (Aqqed), the screened QED effects (A^fsQEo), and nuclear 
effects (A^fnuc)- We thus write the total theoretical value for the g factor as 

g = gY) + Agint + A^Iqed + A^isqed + ^gnnc , (7) 



with the corrections Agi^t, AgqEu, A(7sqed5 and Ag^uc evaluated in Sections |nil|IVl El andlVIl 
respectively. 

III. INTERELECTRONIC INTERACTION 

To evaluate the interelectronic-interaction effects, we start with the relativistic Hamiltonian in 
the no-pair approximation, 

H = A+(J2 ^0) + ^"'') = H^ + H^, (8) 

j 

where h is the one-particle Dirac Hamiltonian the index j = 1 . . . numerates the electrons, 
H'^ is the Coulomb interaction, is the Breit interaction, and A+ is the projector on the positive- 
energy states, which is the product of the one-electron projectors 

A+ = A+(l)---A+(Ar). (9) 

Here, 

A+(z) = 5^|M„(0)(n„(z)|, (10) 

n 

where Un{i) are the positive-energy eigenstates of an effective one-particle Hamiltonian h^. Two 
different Hamiltonians h and h™^ were considered as the operator h^, where h is the one-electron 
Dirac Hamiltonian ^ and h^^^ is the Hartree-Fock-Dirac operator without the Breit interaction 
but including the interaction with the external field Is^Q]- However, since h^^^ corresponds to 
a better zero approximation and, therefore, provides much faster convergence, all the final results 
were obtained with this Hamiltonian. In both cases the functions m„ and hence the projector A+ 
depend on the external magnetic field. Therefore, the derivative with respect to H in Eq. @ con- 
tains not only the derivative of the one-electron part of the Hamiltonian (HI), but also the derivative 
of the projector A+. The derivative of A+ can be obtained explicitly and one can show that its 
contribution to the g factor is equivalent to the contribution of the negative-energy states. 
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In our computational approach, we obtain the (7-factor value by using Eq. ^ in the finite- 
difference approximation, 

since the Hellman-Feynman theorem is not exactly fulfilled for the approximate wave functions. 
The optimal choice of the finite difference was found to be 0.001 (a.u.//io). We checked that 
scaling this value by a factor of 2 does not influence our results for the g factor. 

In order to determine the space of one-electron functions {^Pj}jLi, we employed the combined 
Dirac-Fock (DF) (j = 1, . . . , m) and the Dirac-Fock-Sturm (DFS) (j = m + 1, . . . , M) basis 
set. The index j here enumerates different occupied and vacant one-electron states. The external 
magnetic field and the Breit interaction were not included in the DF and DFS operators, when 
the basis set was generated. For the occupied atomic shells, the orbitals ipj were obtained by 
the restricted DF method, based on a numerical solution of the radial DF equations. The vacant 
orbitals (j' = m + 1, . . . , M) were obtained by solving the Dirac-Fock-Sturm equations 

- e,,] (p^ = \jW{r)(p, , J = m + 1, . . . , M, (12) 

where h^^ is the Dirac-Fock operator, ^j,, is the one-electron energy of the occupied DF orbital 
and W{r) is a constant-sign weight function. The parameter \j in Eq. (fT2b can be considered 
as an eigenvalue of the Sturmian operator. If W{r) ^ at r ^ 00, all Sturmian functions ^pj 
have the same asymptotics at r — > 00. It is clear that for \j = the Sturmian function coincides 
with the reference DF orbital (pj^. The widely used choice of the weight function is W{r) = 1/r, 
which leads to the well-known "charge quantization". In the relativistic case this choice is not very 
successful. In our calculations we used the following weight function 

mr) = ' - 7'Ti°''''' ■ (13) 

which, unlike 1 /r, is regular at the origin. It is well known that the Sturmian operator is Hermitian 
and, contrary to the Fock operator, does not have the continuum spectra. Therefore, the set of 
the Sturmian eigenfunctions including the negative-energy states forms the discrete and complete 
basis set in the space of one-electron wave functions. This basis set is orthogonal with the weight 
function W{r). 

To generate the one-electron wave functions 'tpn, we used the unrestricted DF method in the 
joined DF and DFS basis, 

V'n = ^ Cjn^j . (14) 



The coefficients Cjn were obtained by solving HFD matrix equations 

FCn = EnSCn, (15) 

where F is the Dirac-Fock matrix in the joined basis of DP and DPS orbitals of a free ion. The 
external magnetic field was included in the F matrix, whereas the Breit interaction was not. The 
matrix S in Pq. (fTSb is nonorthogonal, since the DPS orbitals are not orthogonal in the usual 
sense. The negative-energy DPS functions were included in the total basis set. Pq. dTSt was used 
to generate the whole set of orthogonal one-electron wave functions i^nin = 1, • • • M), including 
all vacant states. 

It should be noted that even if the external magnetic field H is equal to zero, the set of one- 
electron functions ipn differs from the set of basis functions (^j. Por the occupied states, the unre- 
stricted DP method accounts for the core-polarization effects (the spin polarization in our case), in 
contrast to the restricted DP method. Por the vacant states the difference is more significant, since 
the DP and DPS operators are essentially different. 

The large-scale configuration-interaction Dirac-Pock (CTDP) method was used to solve the 
Dirac-Coulomb-Breit equation in the external magnetic field 

Hm{-^Mj) = E{-fMj) ^(7Mj) (16) 

where H is the non-pair Hamiltonian The many-electron wave function "^{'jMj) with quan- 
tum numbers 7 and Mj was expanded in terms of a large number of the Slater determinants (SD) 
with the same projection Mj of the total angular momentum J 

^(7Mj) = J2^-{lMj)detUMj). (17) 

a 

The configuration state functions (CSPs) with angular momentum J were not used in our calcu- 
lations, since the Hamiltonian H contains the interaction with the external magnetic field and, 
therefore, does not commute with the operator .P. The Slater determinants are constructed from 
the one-electron wave functions tpn <fT4l) . The same orbitals were used in Pq. (flOb in order to 
construct the projector A+. The basis of one-electron functions used in our calculations was 
12s lip lOd 6/ Ag 2h li. The set of the SD in expansion (flTb was generated including all single, 
double, and triple excitations. The total number of SD was 552359. The results of the calculation 
are presented in Table HI The interelectronic-interaction correction Ag^^~^^ is the difference of 
the result obtained by Pq. dTTT) for the point nuclear model and the Dirac (7-factor value ©. 
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The result for the Ag^^~^^ correction obtained by the CI-DF method can be improved by em- 
ploying a rigorous OED treatment of the part of this correction that is of order 1/Z, which was 
presented in Ref. ll25ll . In order to combine two different treatments, we isolate the contribution 
of order l/Z"^ and higher from the Ag^^^^^ correction by subtracting the value of the 1/Z term 

(2+) 

calculated in the Breit approximation. The resulting "higher-order" correction (A^fj^^^ ) is listed 
in the third column of Table IJ The numerical results for the interelectronic-interaction correction 



of first order in 1 /Z (Agl^l) are taken from Ref. Il25ll and listed in the fourth column. This contri 



bution was evaluated in framework of QED and utilizing the Fermi model for the nuclear-charge 
distribution. Total results A(7int = Ag^l + Ag[l^^ are presented in the last column. The error bars 
indicated represent a quadratical sum of a numerical error and an estimation of omitted terms, i.e., 
contributions beyond the Breit approximation to Agll^\ They were estimated as (aZ)^ Agf^^ 



To compare our results with the corresponding Yan's calculations 10,0]^ which account for 
the lowest-order (~ a^) relativistic effects, we have isolated the contribution in our CI-DF 
calculation. It was done by four times increase of the velocity of light (in atomic units) and by 
an extrapolation of the obtained results (with the factor isolated) to the c = oo limit. Table HH 
which presents the related comparison iox g — 2 values, shows that for Z > 5 the contribution of 
the higher-order relativistic effects is much larger than the difference between our results and 
those of Yan. 



int 



IV. ONE-ELECTRON QED CORRECTIONS 

To zeroth order in the QED correction for the ground state of a Li-like ion is given by the 
one-electron QED contribution evaluated for the 2s Dirac state. This correction is represented by 
a perturbation expansion in the fine-structure constant a, 

A<7QED = A(7«D + Hed + • • • , (18) 

where the superscript indicates the order in a. 

The first-order QED correction is given by the sum of the self-energy and vacuum-polarization 
contributions, A^fq^j-, = Ag^^ + Ag^. The self-energy correction for the 2s state was recently 
calculated to all orders in aZ in Ref. fS"]. The corresponding results are listed in the second column 
of Table Unl The vacuum-polarization correction consists of two parts that can be thought to orig- 
inate from the first-order vacuum-polarization diagram with the magnetic interaction inserted into 
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the external electron line (the electric-loop contribution A^fyp), and into the vacuum-polarization 
loop (the magnetic-loop contribution A^fyp). The electric-loop contribution is calculated in the 
present work by utilizing the known expression for the Uehling potential and approximate for- 



mulas for the Wichmann-KroU potential taken from Ref. 113511 . The results of the calculation are 
presented in the third and fourth columns of Table |ni| The remaining magnetic-loop correction 
is known to vanish in the Uehling approximation, and its contribution is small as compared to the 
electric-loop part. This correction was calculated to all orders in aZ only for the Is state [0,0]. 
Because of this, we employ the analytical result of Ref. Illi for its leading contribution in aZ, 
which reads (for an ns state), 

A«-, = -I-^^M)!. (19) 

216 ri'^ 

The corresponding results are listed in the fifth column of Table |ni| The uncertainty of the 
magnetic-loop contribution is estimated by comparison of the all-order numerical results of 
Ref. lO] for the Is state with the lowest-order analytical result (fT9l) . 

QED corrections of higher orders in a have not been calculated to all orders in a Z up to now. 
Two first terms of their aZ expansion can be re pres ented by the free-electron g—2 factor multiplied 



an can be re pres 



by a relativistic kinematical factor y/ZlliSllZSQa]- The result yields (for an ns state) 

A9« =2n'A«fl + M)!V (20, 



where 2(a/7r)M(*) is the contribution of order to the free-electron g factor. Its numerical values 



(see and references therein) are 



A« = \ , (21) 

^(2) = -0.328 478 965... , (22) 

A^^^ = 1.181241456. . . , (23) 

A(^) = -1.7366(384). (24) 



A.gjh.o.) = {^y^^ {f ln[(aZ)-1 + a.o} , (25) 



Recently, a part of the two-loop QED contribution of order a'^{aZy was evaluated in Ref. 
with the result 



where the numerical values of the coefficient term are ay(ls) = —18.477948 664(1) and 



040^ (2s) = —19.781820 939(1). This expression accounts for the complete logarithmic depen- 
dence in this order and the dominant part of the constant term. We observe, however, that in the 



one-loop case, inclusion of the term of the order [olZY is meaningful for sufficiently small values 
of Z only. In the high-Z region, addition of this term makes the aZ-expansion results deviate 
more from the "exact" numerical values. We thus include the contribution for Z < 30 only. 
The relative uncertainty of the two-loop binding QED correction for Z < 30 was estimated as the 
ratio of the part of the one-loop QED correction that is of order {aZY and higher to the part that 
is within the (aZ)'^ approximation \ multiplied by a factor of 4. For Z > 30, we estimate the 
relative uncertainty by the ratio of the part of the one-loop QED correction that is of order (aZ)^ 
and higher to the part that is within the (aZy approximation, multiplied by a factor of 2. 

V. SCREENED QED CORRECTIONS 

In this section we investigate the influence of the interelectronic interaction on the QED effects, 
known also as the "screening" of QED corrections. To derive the screened QED correction of first 
order in 1/Z, we apply an approximate method, which yields the complete result to the order 
{aZy. Following Hegstrom [39], we adopt the Hamiltonian 

H = J2 Hj) + + ^^^^ 

3 j 

where h is defined by Eq. (O, if™* incorporates the interelectronic interaction within the Breit 
approximation and h^'^'^ accounts for the interaction of the anomalous magnetic moment of 
electron with the magnetic and electric fields, 

h''" = /io [/3(S ■ H) - zP{<x . £)] . (27) 

Here 

^free = 2 5^(^)^A« (28) 



, cr , 

is the free-electron g factor and S = . The magnetic field 7i(j) acting on the jth 

a- 

electron in Eq. (l27t includes the external homogeneous magnetic field Ti and the field induced by 
the other electrons 

n{j) = n + J2T-^^^ (29) 
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where rjk = rj — r^. The electric field £(j) includes the fields induced by nucleus and by the 
other electrons 

^, , \e\Z Tj \ e Vjh 

£U) = ^ 1 + ^1 — • (30) 

We divide the contribution arising from h^'^^ into three parts, 

j 

o-rad _ 9hee - 2 |e|Z Ctj ■ 

^2 - ^ l^Pj—^, (32) 

i/g - /io — 2^ /y^S, ■ —3 ^Pi ^3. I • (33) 



jk jk 



The matrices ctj, Pj, here act on the spinor variables of the jth electron. To first order in 
1/Z, the screened QED correction to the g factor can be now derived by the standard Rayleigh- 
Schrodinger perturbation theory, separating corrections linear in the magnetic field H, which are 
of first order in the parameter 1/Z. The contributions of interest can be conventionally represented 
by the following combinations: 

AEi ~ Hl^'^ X H'''^ + (permutations) , (34) 
AE2 ~ if^^^ X V'^^^'' X H'''^ + (permutations) , (35) 

AE3 ~ Hl'''^ X V"^^"" + (permutations) . (36) 

After angular integration, the summation over the complete Dirac-Coulomb spectrum was per- 
formed by the finite basis set method with basis functions constructed from B splines 14(1 
The numerical results for the screened QED correction of first order in 1/Z, A^fgQp^j-,, are pre- 
sented in Table|IVlin terms of the function R{aZ), defined as 



A^^Sed = (^7frce - 2) ^ RiaZ) . (37) 



The terms Ri{aZ) with i = 1,2, and 3 are induced by Eqs. (I34l) - (l36t . respectively. 

The screened QED correction of higher orders in 1/Z should be accounted for when consider- 
ing low-Z ions. We extract this correction from the recent evaluations by Yan 



30|,|31|], which were 
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performed on nonrelativistic wave functions but with the interelectronic interaction taken into ac- 
count to all orders inl/Z. Yan's results for the screened QED correction can be represented in the 
form 



A5fSQED(Yan) = {gf^^e - 2) {aZf 



(38) 



where -R(O) = —274/2187. The functions R,Q, . . . here do not have any dependence on Z since 
Yan's calculations are based on the nonrelativistic form of the Hamiltonian (l26b . We obtain the 
numerical value of Q{0) by fitting Yan's results for Z = 3-12 to the form (l38t . which yields 
Q{0) = 0.071(1). With this value of Q{0), formula (l38l) was used to estimate the higher-order 
screened QED correction for Z > 12. 

In Table IV] we present the results for the screened QED correction of first order in 1/Z 
(AfyfgQgj-,) and of higher orders in 1/Z (A^fg^j^). The error of the term AgfgQgj) was estimated 
as the part of the one-electron QED correction that is beyond the (aZy approximation, multi- 
plied by a factor 3/Z. The uncertainty ascribed to the A^fg^^^ contribution was evaluated as 
SiaZfAg^^^i^. 

VI. NUCLEAR EFFECTS 

In this section we briefly summarize the known results for the nuclear effects on the g factor 
of Li-like ions. The correction to the Dirac (^-factor value due to the extended nuclear size is 
relatively simple. For high-Z ions, it is evaluated numerically by employing the Fermi model 
for the nuclear-charge distribution. The uncertainty is estimated by taking the difference of the 
results obtained for the Fermi and sphere nuclear models. For \ow-Z ions, to a good accuracy, this 



correction can be evaluated by a simple analytical formula obtained in Ref. lll2ll . 

A,.s - \{aZYm^r^) [l + ^aZf g - C - (:!M^) 



(39) 



Here C = 0.5772156649 ... is the Euler constant and the expectation value has to be evaluated 
with the proper nuclear-charge density. 

Systematic QED theory for the nuclear-recoil effect on the atomic g factor to the first order 
mm/M and to all orders in aZ was developed in Ref. Iil3l1 . The one-electron recoil correction 
derived in that work is expressed as the sum of the lower-order and the higher-order term. The first 
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one can be calculated analytically to yield for the 2s state (cf. 111311 ') 

m{aZf{ 14 + 67 + 12v/2(rT^ , 

^ 4 I ^ ^ 3(1 + ^)2 [2 + ^/2(1+^]' ' 

where 7 = a/1 — (aZ)^. The higher-order term AfyfrecH was calculated numerically in Ref. Iisl 
for the Is state only. These numerical results showed the {aZY behaviour of this term at low Z. 
We estimate the relative uncertainty of the result ( BUt due to neglecting the higher-order term as 
the ratio A^frocH/A^'roc.L for the Is state, multiplied by a factor of 1.5. 

The contribution of the two-electron recoil effect can be extracted from the results of Yan 
HQ. ForZ = 3 - 12, it was evaluated as the difference of Yan's result of order a^m/M and 
the nonrelativistic limit of Eq. (l40b . Fitting this difference to the form 



^^two-el 



'rec 



m {aZf 
M 4 



(41) 



yields C = —3.3(2). With this value of C, formula (1411) was used to estimate the two-electron 
recoil correction for Z > 12. 

finally, we note that the nuclear polarization effect on the atomic g factor was evaluated in Ref. 



VII. RESULTS AND DISCUSSION 

In Table |VT1 we present the individual contributions to the g factor of the ground state of Li- 
like ions. The Dirac point-nucleus value is obtained by Eq. Q. The finite nuclear size correction 
to this value is evaluated by Eq. (l39t for low-Z ions and by a direct numerical solution of the 
Dirac equation for high-Z ions. The interelectronic-interaction correction A^fjnt is the sum of the 
part of first order in 1/Z, AgUl, obtained in framework of QED and of the higher-order part, 
^9iDt^' evaluated by the CI-DF method, as described in Section Unl The QED correction of 
order a is the sum of the one-electron self-energy and vacuum-polarization terms presented in 
Table Uni The QED correction of order and higher incorporates the known terms of the aZ 
expansion, as explained in Section |IV| The screened QED correction is discussed in Section IVl 
the corresponding results are taken from Table [Vj The nuclear recoil correction is obtained as 
explained in Secti on |yil For lead and uranium, we include also the nuclear-polarization correction 
calculated in Ref. M. 
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Table IVTl demonstrates a significant improvement achieved comparing to our previous evalu- 
ations The uncertainty of the presented theoretical values for carbon and oxygen is 2 
and 3 times better than those of Ref. 126], respectively; whereas for uranium the accuracy is im- 
proved by two orders of magnitude. Progress in the high-Z region is mainly due to the evaluation 
of the one-loop QED corrections to all order in aZ, while for \ow-Z ions it is largely due to the 
interelectronic-interaction correction and the screened QED correction. 

The accuracy of the theoretical values in Table I VII is several parts in 10~^ for \ow-Z ions and a 
few parts in 10^^ for middle- Z ions. It decreases further with Z increasing and reaches 5 x 10^^ for 
uranium. So, despite the achieved improvement, the accuracy for Li-like ions is still significantly 
lower than that for H-like ions iQ-Ilsl and also than the precision that can be presently addressed 
in experiments | HO, 3]. In particular, the nearest aim of experimental investigations of the Mainz- 
GSI collaboration is the g factor of Li-like calcium. The anticipated experimental accuracy is at 
the ppb level, which can be compared with the relative theoretical error of 10^^ from Table IVII 

The uncertainty of the present theoretical values is mainly defined by the interelectronic- 
interaction correction and by the screened QED correction. An improvement in the theoretical 
description of the interelectronic-interaction effects can be achieved by a rigorous QED treatment 
of the part of order l/Z"^ and by calculating the remainder within the Breit approximation. Such 

nnnn 

a program has been carried out for the Lamb shift in Li-like ions II42L I43L 1441 14511 . However, a 
similar calculation for the g factor is going to be significantly more difficult due to the presence of 
the external magnetic interaction and requires a further development of methods of calculational 
QED. As to the screened QED correction, in the present work it was calculated to its leading order 
in aZ only. As a first step beyond this approximation, which can improve the results for this cor- 
rection in the high-Z region, one may consider an evaluation of the one-loop QED corrections in 
an effective potential that partly accounts for the interelectronic-interaction effects l^^Hl- These 
two topics will be the subjects of our subsequent investigations. 

In summary, we have presented calculations of the interelectronic-interaction, one-electron 
QED, and screened QED corrections to the g factor of the ground state of Li-like ions. This re- 
sulted in a significant improvement of theoretical predictions in a wide range of the nuclear-charge 
values Z. We have analyzed also perspectives for further progress in the theoretical description of 
these systems and for probing the QED effects in future experiments. 
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The contributions A(/-^^ 
- for the Fermi model of the 
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TABLE 1: The interelectronic-interaction correction A^iint = Ag\Jf, + ^9hit 
and A^;^^"^^ are calculated for the point nucleus, whereas the correction Ag\ 
nuclear-charge distribution. All values shown in units of 10~^. 
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TABLE II: The comparison of the present CI-DF calculations of the g factor and the values obtained by Yan 



|30|,|31j]. The second column presents our g — 2 values, incorporating the effects of binding and electron 
correlation. The limit is shown in the third column. Yan's contribution to the g — 2 value of order is 
listed in third column. All numbers are in units of 10^^ 
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-1012.504 


-1010.072 


-1010.096 
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TABLE III: One-electron self-energy and vacuum-polarization corrections of first order in a. Ai^yp and 
A^'y^^ are the Uehling and the Wichmann-KroU part of the electric-loop vacuum-polarization contribu- 
tion, respectively; A^yp = Ag^yp -|- Ag-y^^ -|- Agiyp. All contributions are calculated with the point 
nuclear model. Units are 10~^. 
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TABLE IV: Individual contributions to the screened QED correction of order 1/Z in terms of the function 
R{aZ) defined by Eq. (ITTl . The terms Ri with i = 1,2,3 are induced by Eqs. (H^-d^^. respectively. 
The results in the last column are obtained for the extended-charge nucleus, whereas all the other results 
correspond to the point nuclear model. 
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TABLE V: The total screened QED correction. A(7gQgj^ and Ag'gQ^j-, are the screened QED contribution 
of first and of lugher orders in 1/Z, respectively. All numbers are in units of 10~^. 
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TABLE VI: Individual contributions to the ground-state g factor of Li-Uke ions. 
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